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Abstract 

The  main  objective  of  this  work  is  to  extract  information  on  the  internal  behavior  of  three  small  polymer  electrolyte  membrane  fuel  cells  under 
static  and  dynamic  load  conditions.  A  computational  model  was  developed  using  Scilab  [SCILAB  4,  Scilab — a  free  scientific  software  package, 
http://www.scilab.org/,  INRIA,  France,  December,  2005]  to  simulate  the  static  and  dynamic  performance  [J.M.  Correa,  A.F.  Farret,  L.N.  Canha, 
An  analysis  of  the  dynamic  performance  of  proton  exchange  membrane  fuel  cells  using  an  electrochemical  model,  in:  27th  Annual  Conference 
of  IEEE  Industrial  Electronics  Society,  2001,  pp.  141-146]  of  this  particular  type  of  fuel  cell.  This  dynamic  model  is  based  on  electrochemical 
equations  and  takes  into  consideration  most  of  the  chemical  and  physical  characteristics  of  the  device  in  order  to  generate  electric  power.  The  model 
takes  into  consideration  the  operating,  design  parameters  and  physical  material  properties.  The  results  show  the  internal  losses  and  concentration 
effects  behavior,  which  are  of  interest  for  power  engineers  and  researchers. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

In  the  last  decade,  a  group  of  researchers  [5,6]  published  a 
series  of  papers  on  electrochemical  problems  relating  to  PEM 
fuel  cells,  from  modeling  to  the  problems  of  poisoning  of 
the  cell,  their  work  aimed  at  the  need  for  a  model  based  on 
electrochemical  equations  to  correctly  model  this  new  elec¬ 
trical  power  source.  Almost  10  years  later  another  group  of 
researchers  has  focused  on  electrochemical  reaction  model¬ 
ing  [4],  but  like  previous  publications,  their  model  is  based 
on  the  assembling  of  a  prototype  fuel  cell  and  laboratory 
experiments. 

More  recent  work  by  Souza  and  Gonzales  [15]  has  presented 
a  mathematical  model  for  the  PEM  Fuel  Cell,  and  their  modeling 
has  allowed  detailed  studies.  They  discussed  the  electrocatalysis 
of  the  reactions  and  the  design  of  water-management  schemes  to 
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cope  with  membrane  dehydration.  Mishra,  Yang  and  Pitchumani 
[18],  did  a  detailed  numerical  investigation  of  the  transport  and 
electrochemical  phenomena  involved  in  the  operation  of  a  sin¬ 
gle  proton  exchange  membrane  (PEM)  fuel  cell  with  reformate 
feeding  with  a  view  to  developing  optimal  design  and  operating 
conditions. 

Along  with  a  large  number  of  methodologies,  Cheddie  and 
Munroe  [16]  did  a  review  of  recent  literature  on  proton  exchange 
membrane  fuel  cell  modeling,  and  they  have  categorized  them 
as  analytical,  semi-empirical  or  mechanistic.  Mechanistic  mod¬ 
eling  has  received  the  most  attention  in  the  literature.  In  mecha¬ 
nistic  modeling,  differential  and  algebraic  equations  are  derived 
based  on  the  physics  and  electrochemistry  governing  the  phe¬ 
nomena  internal  to  the  cell. 

Mehta  and  Cooper  [17]  did  an  extensive  review  and  analysis 
of  PEM  fuel  cell  designs  and  manufacturing.  Their  work  had 
an  important  influence  on  the  work  we  present  here,  Goveli- 
oglu  and  Stenger  [19]  developed  research  on  computational  fluid 
dynamics  modeling  of  PEM  fuel  cells.  They  used  a  finite  ele¬ 
ment  method  to  solve  the  transport  model  coupled  with  flowing 
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Nomenclature 

A 

anode  (or  cathode)  area/cell  useful  area  (cm2) 

* 

ch2 

hydrogen  concentration  in  the  interface  electrode 
membrane  (mol  cm-3) 

Cion,m 

membrane  ionic  concentration  (mol  cm-3) 

£° 

open  circuit  voltage  (V) 

F 

Faraday  constant  (96487  Cmol-1) 

AFEa 

anode  reaction  activation  energy  (J  mol-1) 

AFEC 

cathode  reaction  activation  energy  (J  mol-1) 

AG 

Gibbs  free  energy  variation  (—237.165  J) 

AH 

variation  of  enthalpy  (—285.823  J) 

lop 

operational  cell  current 

ih 

electrode  current  limit 

A, a 

electrode  current  limit  in  the  anode 

*L,c 

electrode  current  limit  in  the  cathode 

lop 

operational  electrical  current  (A) 

7op 

cell  operation  current  density  (A  cm-2) 

jo 

exchange  current  density  in  the  electrode 
(A  cm-2) 

k° 

reaction  degree  intrinsic  constant  (cms-1) 

l 

membrane  thickness  (cm) 

PlJ\r 

mass  flow  of  air  (dimensionless) 

mu2 

mass  flow  of  hydrogen  (dimensionless) 

mo2 

mass  flow  of  oxygen  (dimensionless) 

n 

number  of  moles  involved  in  a  reaction  for  each 
mol  of  H2  (dimensionless) 

nSL 

number  of  moles  of  reagent  in  the  anode 

nc 

number  of  moles  of  reagent  in  the  cathode 

Pr2 

partial  pressure  feeding  for  hydrogen  (atm) 

Po2 

partial  pressure  feeding  for  oxygen  (atm) 

R 

universal  constant  of  gases  (8314  JKmol-1) 

Rcj 

resistance  of  gas  spreaders  and  separating  plates 
(«) 

R\u 

resistance  of  the  electrodes  (Q) 

AS 

variation  of  entropy  (  — 163.2  J  K-1) 

T0p 

operational  temperature  (K) 

Tref 

reference  temperature  (298.15  K-1) 

katv 

activation  voltage  drop  (V) 

Econ 

concentration  voltage  drop  (V) 

Eohm 

Ohmic  voltage  drop  (V) 

Fop 

operational  voltage  (V) 

Greek  letters 

a 

electronic  transference  current  (around  0.4  for  the 
cathode  and  1 .0  for  the  anode) 

ac 

cathode  load  transference  coefficient  (dimension¬ 
less) 

Ota 

anode  load  transference  coefficient  (dimension¬ 
less) 

y 

membrane  ionic  conductance  (S  cm2  mol-1) 

Sc 

cathode  Nernst  diffusion  layer  thickness  (cm) 

5a 

anode  Nernst  diffusion  layer  thickness  (cm) 

^-air 

stoichiometry  of  air  (dimensionless) 

■^h2 

stoichiometry  of  hydrogen  (dimensionless) 

Ao2  stoichiometry  of  oxygen  (dimensionless) 

pm  electrical  resistivity  of  the  membrane  (Q  cm) 
0H?  relative  humidity  of  hydrogen  (%) 

0o2  relative  humidity  of  oxygen  (%) 


in  a  porous  medium,  charge  balance,  electrochemical  kinetics, 
and  a  rigorous  water  balance  in  the  membrane. 

These  papers  have  contributed  greatly  to  comprehension 
of  the  global  electrochemical  processes  in  a  fuel  cell,  but 
did  not  improve  understanding  of  the  fuel  cell  as  an  electri¬ 
cal  energy  source,  and  that  is  our  proposal  in  this  work,  i.e. 
to  make  a  model  [1]  suitable  for  use  with  power  electronics 
hardware. 

The  basic  physical  structure  of  a  fuel  cell  used  in  the  computer 
model  consists  of  an  electrolyte  and  a  polymeric  membrane,  with 
a  typical  thickness  of  between  12  and  210  p,m,  which  separates 
the  porous  anode  and  cathode  from  each  other;  they  have  a  thick¬ 
ness  of  between  5  and  15  pan.  The  porous  and  slender  electrodes, 
make  it  possible  to  assemble  a  membrane  electrode  assembly 
(ME A).  In  a  typical  fuel  cell,  a  flux  of  fuel  which  is  rich  in  hydro¬ 
gen  (liquid  or  gaseous)  is  fed  continuously  through  channels  that 
keep  it  in  contact  with  the  anode  and,  simultaneously,  a  flux  of 
an  oxidant  (in  general  the  oxygen  in  the  air)  are  also  fed  through 
different  and  separate  channels,  which  make  the  contact  with  the 
cathode  [2].  The  layers  in  between  are  responsible  for  the  gas 
diffusion  process  in  the  electrodes,  with  a  thickness  between  300 
and  400  [xm,  they  are  made  of  an  electrical  conducting  material, 
like  porous  carbon  paper,  a  carbon  based  fabric  with  a  Teflon® 
layer  that  prevents  water  adherence  and  make  possible  the  fast 
diffusion  of  gases  [3]. 

The  process  that  takes  place  is  a  natural  reaction,  however 
due  to  its  slow  speed  it  is  necessary  to  use  a  catalyst  for  the 
reaction  to  be  efficient,  the  element  used  as  catalyst  at  the  present 
is  platinum.  The  other  two  important  elements  in  the  fuel  cell 
are,  firstly  the  separating  plates,  to  give  mechanical  stiffness 
and  to  direct  the  flow  of  gas  through  the  cell,  and  secondly  the 
cooling  elements  that  have  the  function  of  extracting  heat  from 
the  exothermal  reactions,  the  ratio  used  for  interspersing  the 
cooling  elements  are  used  in  a  proportion  that  varies  from  1 : 1 
to  1:5  [5]. 

2.  Mathematical  modeling  of  the  PEM  fuel  cell 

A  mathematical  model  that  presents  the  correct  answers  for 
the  fuel  cell  is  of  paramount  importance,  since  it  allows  the  devel¬ 
opment  of  electronic  systems  for  feeding  and  control  of  fuel  cells 
without  the  need  for  a  great  number  of  experimental  prototypes. 
Electrically  speaking,  a  fuel  cell  works  basically  as  a  current 
source.  The  voltage  response  for  the  cell  will  be  studied  from 
the  current  variation  required  by  the  load.  The  main  operational 
parameters  are:  the  operational  electrical  current,  temperature, 
and  the  characteristics  of  gas  transportation. 
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2.1.  Activation  losses 


2.2.  Resistive  losses 


2.1.1.  Cathode  activation  losses 

The  consideration  on  the  effect  of  operating  variables  in  the 
value  of  io  in  the  cathode  can  be  taken  as  [6] : 


ioc  =  ncFAk\C*pmtof 


-ac 


Working  mathematically  in  the  previous  equation  we  can  have 
for  the  cathode  the  following  value  for  the  activation  voltage  [4] : 


fatvc  —  +  $2Top  +  &  Tip  ln(Cj2)  +  $4T0 p  ln(i) 


(2) 


where 


Si  = 


$2  = 


S3  = 


&  = 


AFEC 

nca  c 

R 

ncFac 


ln(4  Fk°)  +  — 7 — ln(A)  +  ^ln(C^) 


R 


R 


ncFac 


ricF 


Rjl  -ac) 
ncFac 

R 


ncFac 


(3) 


(4) 


(5) 


(6) 


2.1.2.  Anode  activation  losses 

With  a  more  detailed  Eq.  ( 1 ),  taking  in  consideration  the  effect 
of  operational  variables  in  the  value  of  io  at  the  anode  as  [6] : 

io,a  (A)  =  n3FAk°Cllf-AFE^RT^>  (7) 

The  activation  voltage  [4]  in  the  cathode  can  be  given  by: 


Vatv,c  (V)  =  §5  +  ?6^op  ln(C*2)  +  ^op  ln(0 


where 


§5  = 


$6  = 


A  FE 


a 


naFaa 


-RXnfn^FAk0) 


nzFcta 


Hi  = 


R 


n^Foi^ 


(8) 

(9) 

(10) 

(11) 


As  previously  stated,  the  activation  voltage  at  the  anode  will 
have  lower  values  compared  with  the  activation  voltage  at  the 
cathode,  and  can  be  neglected  in  the  operating  behavior  of  the 
PEM  fuel  cell  fed  with  pure  hydrogen. 


To  illustrate  the  empirics  to  obtain  fuel  cell  operating  param¬ 
eters,  the  first  part  of  the  resistive  losses  (R0hm)  Is  given  by 
the  polymeric  membrane  and  the  second  part  by  the  electrodes, 
spreaders  and  other  series  components. 


2.2.1.  Polymeric  membrane  resistive  losses 

The  resistance  in  Ohm  for  this  membrane  ( Rm )  which  rep¬ 
resents  the  resistance  for  transference  of  protons  through  the 
solid  membrane  [14],  and  can  be  computed  by  the  two  follow¬ 
ing  equations: 

Pm  l 

Rm(£2)=^  (12) 

The  resistivity  values  for  the  Nation®  [9]  membranes  (produced 
by  DuPont)  to  complement  the  ohmic  resistance  can  be  com¬ 
puted  through  the  following  equation: 


pm  (Q  cm)  = 


181.6 

1  +  0.03(iop/A)  +  0.0062(Top/303)2(7op 

/a)2-5; 

A  —  0.634  —  3(/Op/A)e4-18((7op_303)/(:rop)) 

(13) 


For  the  Nation®  membrane,  the  typical  thickness  in  each  case 
is  indicated  in  Table  1  [8,12]. 

Therefore,  X  is  the  ratio  of  the  number  of  water  moles  for 
each  sulfonic  group  in  the  membrane.  The  value  of  A  indicates 
the  membrane  relative  humidity,  being  regulated  by  the  feeding 
gases  relative  humidity,  their  stoichiometries  and  assembling 
characteristics  involved  in  the  efficiency  of  the  gas  transportation 
and  membrane  humidity  maintenance.  The  literature  indicates  a 
value  of  approximately  14  for  A  corresponding  to  a  membrane 
in  humidified  air  at  100%,  and  approximately  to  22  the  value 
corresponding  to  a  water  immersed  membrane  [4] . 


2.3.  Concentration  losses 


2.3.1.  Concentration  losses  at  the  cathode 

The  value  for  the  current  limit  at  the  cathode  (/l,c)  will  be 
adopted  [13]  through  the  following  equation: 


iUc  (A)  = 


RT0pAyC{ 


ion,m 


nrF8( 


(14) 


2.3.2.  Concentration  losses  at  the  anode 

They  will  be  simulated  as  in  Eq.  (14),  and  the  value  for  the 
current  limit  at  the  anode  (/L,a)  will  be  adopted  through  the  fol- 


Table  1 

Typical  thickness  of  Nafion®  membranes 


Membrane 

Nominal  thickness  (mm) 

Thickness  when  dry  (mm) 

Thickness  when  wet  (mm) 

Thickness  when  ME  A  (mm) 

Nafion  117 

0.178 

0.183  ±  0.003 

0.208  ±  0.005 

0.148  ±  0.002 

Nafion  115 

0.127 

0.141  ±  0.003 

0.161  ±  0.003 

0.100  ±  0.002 

Nafion  113 

0.089 

0.091  ±  0.002 

0.111  ±  0.002 

0.075  ±  0.003 

Nafion  112 

0.051 

0.050  ±  0.002 

0.058  ±  0.003 

0.040  ±  0.002 
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0 


temperature  meter 


humidifier 


Fig.  1.  Typical  elemental  operating  scheme  for  a  PEM  fuel  cell  stack. 


lowing  equation: 

R  Tov  A  yCion?m 


A, a  (A)  = 


^a^a 


2.4.  PEMFC  cell  under  dynamic  load  change 


(15) 


When  a  variation  occurs  in  the  load  fed  by  a  fuel  cell,  this  is 
the  cause  for  a  variation  in  the  cell  operating  voltage.  However, 
a  transient  state  exists  due  to  the  capacitance  generated  in  the 
polymeric  membrane  due  to  an  internal  physical  process  called 
“double  layer  effect”.  [7] 

The  equation  that  represents  the  voltage  behavior  during  load 
transients  can  be  expressed  as: 

v(t)(V)  =  v(oo)  +  [p(0+)  —  n(oo)]exp  ^ ^  (16) 

where  v(oo)  is  the  final  value  to  which  the  voltage  will  go  after 
the  load  variation,  and  n(0+)  is  the  voltage  initial  value,  before 
the  variation  of  the  load  fed  by  the  PEM  fuel  cell.  The  time 
constant  (r)  is  a  measure  of  the  necessary  time  for  the  transient 
analysis  in  RC  circuits,  computed  by  the  following  equation: 

r(^)  =  RthCdc,  (17) 

and  Rth  is  the  Thevenin  equivalent  resistance  of  the  circuit  seeing 
from  the  capacitance. 


3.  Methods 


membrane),  and  the  dimensioning  parameters  for  the  fuel  cell 
stack  (number  of  cells,  cell  area  and  membrane  thickness)  were 
also  considered.  The  simulation  results  are  compared  with  prac¬ 
tical  experiments. 

The  experimental  tests  were  made  through  controlled  pro¬ 
cesses  using  different  models  and  prototypes  of  PEMFC  fuel 
cell  stack  as  the  system  shown  in  the  basic  scheme  of  Fig.  1. 

Under  normal  operation  a  frequency  inverter  is  used  to  feed 
ac  energy  due  to  the  voltage  regulation  showed  by  the  fuel  cell. 
An  electronically  controlled  resistive  load  was  used  in  order  to 
make  possible  a  better  control  of  the  necessary  electrical  load 
during  laboratory  experiments. 

4.  Results 

4.1.  PEM  fuel  cell  with  a  load  in  steady-state 

Fig.  3  shows  the  behavior  of  a  cell  with  net  area  of  5  cm2, 
fed  by  hydrogen  and  pure  oxygen  with  feeding  pressures  of 
1  atm,  stoichiometries  of  1 . 1 ,  humidification  100%  under  operat¬ 
ing  temperatures  of 40, 60  and  80  °C  for  a  defined  loading,  being 
the  curves  assembled  from  the  cell  voltage  response  according  to 
the  current  requested  by  the  load.  Similar  curves  can  be  obtained 
fixing  the  temperature  and  changing  other  operating  parameters. 

Fig.  2  shows  the  behavior  of  the  cell  with  the  variation  of 
temperature  for  the  first  prototype  with  constants  described  in 
Table  2  [20]: 


In  order  to  develop  this  study,  computational  modeling  tools 
were  built  up  to  model  the  fuel  cell  behavior  through  a  detailed 
theoretical  base  taking  into  consideration  assembling  and  oper¬ 
ating  parameters. 

With  this  purpose;  to  build  the  fuel  cell  simulator  besides 
the  operating  parameters  already  mentioned  (Top ,  7op ,  Pu2 ,  Pq2  , 
0h2>  0o2  ’  ^h2,  ^o2>  or  A. <rir),  the  constructive  parameters  (elec¬ 
tronic  transference  coefficients,  electrode  reaction  level  con¬ 
stants,  and  the  physicochemical  properties  of  electrodes  and 


Table  2 


Design  characteristics  for  the  first  prototype 


Description 

Symbol 

Amount 

Area  of  the  cell 

A  cel 

5  cm2 

Number  of  cells 

A%1 

1 

Membrane 

— 

Nafion®  112 

Electronic  transfer  coefficient  (cathode) 

ac 

0.43 

Electronic  transfer  coefficient  (anode) 

aa 

1.0 

Reaction  grade  constant  (cathode) 

*? 

2.7  cms-1 
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Fig.  2.  Behavior  of  a  PEM  fuel  cell  under  static  load  with  the  temperature. 


Fig.  3.  Behavior  of  a  PEM  fuel  cell  under  static  load  with  the  gas  pressure. 


Fig.  3  shows  the  behavior  of  the  cell  with  the  variation  of 
pressure  for  the  second  prototype  with  constants  described  in 
Table  3  [20]: 

The  third  prototype  was  used  for  experimental  and  dynamic 
simulation  tests,  it  is  a  fuel  cell  stack  built  from  two  cells,  assem¬ 
bled  in  a  test  bench  with  main  design  characteristics  as  shown 
Table  4  [20]. 


Table  3 


Design  characteristics  for  the  second  prototype 


Description 

Symbol 

Amount 

Area  of  the  cell 

ACel 

20  cm2 

Number  of  cells 

/Vccl 

1 

Membrane 

— 

Nafion®  112 

Electronic  transfer  coefficient  (cathode) 

ac 

0.43 

Electronic  transfer  coefficient  (anode) 

Oi  a 

1.0 

Reaction  grade  constant  (cathode) 

*c° 

3.5  cms-1 

Table  4 


Design  characteristics  for  the  third  prototype 


Description 

Symbol 

Amount 

Area  of  the  cell 

A  cel 

80  cm2 

Number  of  cells 

A^CXl 

2 

Membrane 

— 

Nafion®  112 

Electronic  transfer  coefficient  (cathode) 

O'  c 

0.43 

Electronic  transfer  coefficient  (anode) 

a  a 

1.0 

Reaction  grade  constant  (cathode) 

1.3  cms-1 

4.2.  Utilization  of  air  instead  of  pure  oxygen 

For  this  analysis  the  second  prototype  was  used  with  air  or 
oxygen  as  a  cathode  reagent  with  an  operating  temperature  of 
50  °C.  The  other  cell  parameters  are  the  ones  specified  in  Table  2. 
The  values  obtained  experimentally  and  through  computational 
simulation  are  indicated  in  Fig.  4  [20]. 

The  result  above  shows  the  developed  model  efficiency  and 
the  success  in  the  collection  of  constructive  parameters,  and  also 
demonstrates  the  influence  of  the  use  of  atmospheric  air  and  pure 
oxygen. 

Practically,  this  substitution,  although  showing  loss  of  effi¬ 
ciency,  is  advantageous  for  commercial  projects,  avoiding  the 
need  for  oxygen  transport  and  storage,  since  one  of  the  major 
challenges  for  the  PEM  fuel  cells  is  to  find  an  efficient  way  to 
store  bulky  amounts  of  gaseous  fuels  like  hydrogen  and  oxygen. 

The  decreasing  efficiency  is  due  to  the  lower  oxygen  con¬ 
centration  in  the  composition  of  atmospheric  air,  and  to  the 
consequent  increase  in  the  cathode  concentration  effect  (Fig.  5) 
[20] .  The  ohmic  voltage  drop  and  activation  effect  are  not  depen¬ 
dent  on  this  parameter. 

A  second  reason  for  this  variation  is  the  concentration  effect 
on  the  cathode.  With  a  decrease  in  concentration  of  the  fuel  gas, 
the  maximum  possible  current  in  this  electrode  decreases  [11]. 
Additionally  the  experiment  was  at  a  relatively  low  temperature 
(50  °C),  which  makes  the  cell  response  curve  for  high  current 


Fig.  4.  Behavior  of  a  PEM  Fuel  Cell  for  pure  oxygen  and  air. 
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Ohmic  and  Activation  Effects 


Fig.  5.  Ohmic  and  activation  losses  when  comparing  feeding  of  pure  oxygen 
and  air. 


Concentration  effect 


densities  close  to  the  saturation  limit  (the  oxygen  concentration 
in  the  membrane-electrode  interface  is  close  to  zero).  This  phe¬ 
nomenon  occurs  in  a  region  where  the  power  is  decreasing. 

Fig.  6  shows  also  that  the  anode  concentration  effect  is  mini¬ 
mal,  and  due  to  this  very  small  value,  is  irrelevant  to  the  process 
[20]. 

5.  Conclusion 

As  can  be  seen  from  this  investigation,  the  coupling  of  a  fuel 
cell  with  an  inverter  to  produce  ac  current  is  not  as  simple  as 
it  seems.  The  fuel  cell  as  an  energy  source  changes  its  internal 
parameters  dynamically  for  each  load  or  operating  point.  The 
inverter  design  is  not  simple  since  its  characteristics  should  be 
in  agreement  with  the  fuel  cell  variable  parameters. 

The  study  of  the  dynamic  behavior  shows  that  the  electrical 
RC  constant  can  change  considerably  with  the  control  setup  for 


the  stack  with  all  internal  parameters  and  inputs,  and  also  shows 
that  for  a  certain  operating  point,  the  RC  constant  during  the 
load  increase  is  larger  than  during  the  load  decrease;  opposite  to 
the  fact  that  the  resistance  in  general  increases  with  the  value  of 
current  and  temperature. 

The  implementation  of  this  model  allows  the  assembling  of 
different  sizes  of  PEMFCs,  and  shows  the  effects  of  variations 
on  the  many  components  of  the  performance  of  the  cells. 

Anode  concentration  and  activation  effects  proved  to  be  much 
less  dependent  of  the  type  of  feeding  (air  or  pure  O2)  than  the 
effects  at  the  cathode,  as  can  be  seen  in  Fig.  6.  Ohmic  volt¬ 
age  drops  have  a  non-linear  dependency  with  current  and  can 
be  smaller  than  the  anodic  activation  losses  for  a  given  operat¬ 
ing  current  density. The  developed  computational  simulator  was 
demonstrated  to  be  an  efficient  and  valuable  tool  in  the  develop¬ 
ment  of  polymeric  membrane  fuel  cell  stacks;  making  possible 
an  estimation  of  the  cell  performance  [10]. 


Acknowledgment 

The  authors  wish  to  thank  to  UNITECH  Ltd.,  represented  by 
Dr.  Antonio  Cesar  Ferreira,  Director  of  Research  and  Develop¬ 
ment. 


References 

[1]  SCILAB  4,  Scilab — a  free  scientific  software  package,  http://www.scilab. 
org,  INRIA.  France,  December,  2005. 

[2]  J.H.  Hirschenhofer,  D.B.  Stauffer,  R.R.  Engleman,  M.G.  Klett,  Fuel 
Cell  Handbook,  5th  ed.,  Parsons  Corporation,  Reading/PA,  for  the  U.S. 
Department  of  Energy  Office  of  Fossil  Energy,  PA,  West  Virginia, 
2000. 

[3]  S.  Thomas,  M.  Zalbowitz,  Fuel  Cells — Green  Power,  Los  Alamos  National 
Laboratory,  U.S.  Department  of  Energy,  Office  of  Advanced  Automotive 
Technologies,  2001. 

[4]  J.M.  Correa,  A.F.  Farret,  L.N.  Canha,  An  analysis  of  the  dynamic  perfor¬ 
mance  of  proton  exchange  membrane  fuel  cells  using  an  electrochemical 
model,  in:  27th  Annual  Conference  of  IEEE  Industrial  Electronics  Society, 
2001,  pp.  141-146. 

[5]  J.C.  Amphlett,  R.F.  Mann,  B.A.  Peppley,  PR.  Roberge,  A.  Rodrigues,  A 
Practical  PEM  Fuel  Cell  Model  for  Simulating  Vehicle  Power  Sources, 
Department  of  Chemistry  and  Chemical  Engineering,  Royal  Military  Col¬ 
lege  of  Canada,  Ontario,  1995,  pp.  221-226. 

[6]  A.  Rodrigues,  J.C.  Amphlett,  R.F.  Mann,  B.A.  Peppley,  PR.  Roberge, 
Carbon  Monoxide  Poisoning  of  Proton-Exchange  Membrane  Fuel  Cell, 
Department  of  Chemistry  and  Chemical  Engineering,  Royal  Military  Col¬ 
lege  of  Canada,  Ontario,  1995,  pp.768-773. 

[7]  A.R.  Balkin,  Modeling  a  500  W  Polymer  Electrolyte  Membrane  Fuel  Cell, 
Faculty  of  Engineering,  University  of  Technology,  Sydney,  Australia,  2002. 

[8]  S.  Slade,  S.A.  Campbell,  T.R.  Ralph,  F.C.  Walsh,  Ionic  conductivity  of  an 
extruded  Nafion  1100  EW  series  of  membranes,  J.  Electrochem.  Soc.  49 
(2002)  1556-1564. 

[9]  T.  Suzuki,  H.  Murata,  T.  Hatanaka,  Y.  Morimoto,  Analysis  of  the  Catalyst 
Layer  of  Polymer  Electrolyte  Fuel  Cells,  R&D  Review  of  Toyota  CRDL, 
volume  39,  number  3,  2003. 

[10]  J.B.  Benziger,  et  al.,  The  Auto  Humidification  Polymer  Electrolyte  Mem¬ 
brane  Fuel  Cell,  Princeton  University,  Princeton,  2001. 

[11]  R.R.  Passos,  E.A.  Ticianelli,  Effects  of  the  operational  conditions  on  the 
membrane  and  electrode  properties  of  a  polymer  electrolyte  fuel  cell,  J. 
Braz.  Chem.  Soc.  4  (2002)  483-489. 


1182 


R.A.  Costa,  J.R.  Camacho  /  Journal  of  Power  Sources  161  (2006)  1176-1182 


[12]  S.  Slade,  et  al.,  Ionic  conductivity  of  an  extruded  Nafion  1100  EW  Series 
of  membranes,  J.  Electrochem.  Soc.  149  (2002)  1556-1564. 

[13]  P.R.  Pathapati,  X.  Xue,  J.  Tang,  A  new  dynamic  model  for  predicting  tran¬ 
sient  phenomena  in  PEM  fuel  cell  systems,  Elsevier  Renewable  Energy, 
2004. 

[14]  P.W.  Atkins,  J.  de  Paula,  Physical  Chemistry,  7th  ed.,  Oxford  University 
Press,  2002. 

[15]  R.  Souza,  E.R.  Gonzales,  Mathematical  modeling  of  polymer  electrolyte 
fuel  cells,  J.  Power  Sources  147  (1-2)  (2005)  32-45. 

[16]  D.  Cheddie,  N.  Munroe,  Review  and  comparison  of  approaches  to  proton 
exchange  membrane  fuel  cell  modelling,  J.  Power  Sources  147  (1-2)  (2005) 
72-84. 


[17]  V.  Mehta,  J.S.  Cooper,  Review  and  analysis  of  PEM  fuel  cell 
design  and  manufacturing,  J.  Power  Sources  114  (1)  (2003)  32- 
53. 

[18]  V.  Mishra,  F.  Yang,  R.  Pitchumani,  Analysis  and  design  of  PEM  fuel  cells, 
J.  Power  Sources  141  (1)  (2005)  47-64. 

[19]  G.H.  Guvelioglu,  H.G.  Stenger,  Computational  fluid  dynamics  of  polymer 
electrolyte  membrane  fuel  cells,  J.  Power  Sources  147  (1-2)  (2005)  10-95. 

[20]  R.A.  Costa, Computational  analysis  and  modeling  of  PEM  fuel  cells  under 
static  and  dynamic  load  regimen,  Electrical  Engineering  MSc  Dissertation, 
Universidade  Federal  de  Uberlandia,  2005,  Uberlandia,  MG,  Brazil  (in 
Portuguese). 


